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INTRODUCTION 


I . 

This  report  addresses  three  problems.  First,  new  results  are 
presented  on  the  problem  of  linear  least-squares  restoration  of  photon- 
limited,  atmospherically  degraded  images.  Second,  the  photon-limited 
performance  of  a system  which  combines  pre-detection  compensation  and 
post-detection  compensation  is  considered.  Thirdly,  we  present  some 
recently  obtained  results  on  the  photon-limited  sensitivity  of  the 
"speckle  imaging"  technique  for  obtaining  high  resolution  object  inform- 
ation in  the  presence  of  atmospheric  turbulence. 

The  analysis  and  results  presented  in  this  report  depend  heavily 
on  material  published  in  a previous  technical  report  [1].  As  far  as 
possible,  the  notation  in  this  report  is  identical  to  that  used  in  the 
former,  except  for  a few  small  improvements.  We  shall  feel  free  to  call 
upon  the  previously  derived  results  as  needed. 

I I . BACKGROUND 

The  problem  addressed  here  is  illustrated  in  Fig  1.  A distant 
object,  represented  by  a radiance  distribution  o(x,y)  , is  subjected  to 


Fig.  1.  Block  diagram  for  least-mean-square-error  restoration 


a linear,  isoplanatic  blur,  represented  by  optical  transfer  function 

(The  spatial  frequencies  £2^  and  fty  are  measured  in  cycles 
per  radian  of  arc  in  the  sky.)  In  the  case  of  special  interest  here, 
this  blur  arises  from  both  atmospheric  inhomogeneities  and  diffraction 
due  to  the  finite  aperture  of  the  telescope  used.  The  blurred  irradiance 
distribution  i (x,y)  is  then  detected.  The  detected  signal  d(x,y)  is 
represented  as  a compound  Poisson  process,  with  space  variant  mean  A(x,y) 
proportional  to  i(x,y).  The  detected  signal  is  passed  through  a linear 
space- i nvar  iant  restoration  filter  with  transfer  function  H(ft^,£2y). 

The  restoration  filter  is  chosen  to  minimize  the  mean-squared 
value  of  the  error  e(x,y)  , which  represents  the  difference  between  the 
restored  image  r(x,y)  and  an  “ideally"  filtered  version  of  the  object 
o(x,y).  The  ideal  filter  is  represented  by  a transfer  function  S(ftx,fiy). 

The  quality  of  the  restored  image  is  represented  by  two  parameters. 
The  “restored  bandwidth",  Aft  , is  defined  by 


An  = 


-l 


BHjdn  dn 


J 
— co 


X y 


(I) 


where  B and  H are  both  normalized  to  have  value  unity  at  the  origin, 
The  mean-squared  signal -to-nolse  ratio  in  the  restored  image  is  defined 
by  [t,  pp. 16,17] 


(N)2 


s|2  4>o  dnxd£2y 


(2) 


where  N represents  the  mean  number  of  photoevents  detected  in  the  image, 
$o  is  the  spectral  density  of  the  object,  and  e is  the  expected  value 
of  the  total  mean-squared  error  in  the  restored  Image.  We  have  previously 
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shown  that  e is  given  explicitly  by  [1,  p . 1 6] 
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The  normalized  spectral  density  of  the  object,  , was  shown 

[1  , pp. 18-21]  to  depend  on  the  complexity  of  the  object.  In  the  case  of 
an  object  representable  as  a windowed  sample  function  of  a stationary 
random  process,  the  normalized  spectral  density  was  found  to  drop  rapidly 
from  unity  to  a level  (aQ/o)  (A  /A)  , where  oQ  and  o are,  respectively, 
the  standard  deviation  and  mean  of  the  object  radiance  distribution,  A 

c 

is  the  correlation  area  of  the  object,  and  A is  the  area  of  the  total 
object  field.  In  the  results  to  be  presented  later,  we  shall  always 
assume  that  the  object  is  a point  source  ($>Q  * 1).  However,  to  a good 
approximation  the  results  can  be  applied  to  more  complicated  objects 
provided  the  parameter  N is  interpreted  as  being  equal  to  (aQ/o)  n , 
where  n represents  the  average  number  of  photoevents  arising  from  a 
single  correlation  area  of  the  object. 


III.  RESTORATION  OF  ATMOSPHERICALLY  DEGRADED  IMAGES 

In  our  previous  report  we  showed  plots  of  AH  and  Q vs.  N 
for  values  of  rQ  (the  coherence  length  of  the  atmospheric  distortions) 
of  5cm  and  10cm,  and  for  a flat  ideal  transfer  function 


s(r<x,ny) 


0 otherwise 


(*> 
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where  S2Q  is  the  d i f f ract ion- 1 imi ted  cutoff  frequency  of  the  telescope. 
We  have  extended  the  calculations  to  include  in  addition  the  cases  of 
rQ  = 20cm  and  rQ  = “ (diffraction  limited),  and  we  have  modified  our 
choice  of  the  ideal  transfer  function,  using  instead  the  diffraction- 
limited  transfer  function, 


2 J -1  f2  Si 

S Si)  = - cos  — - — 

* Si  Si 


(2,  2 \ 

+ qy)  2 ’ 

Our  reason  for  changing  the  choice  of  S lies  in  a pathological 
behavior  associated  with  the  flat-topped  S.  When  N has  reached  such  a 
large  value  that  the  restored  bandwidth  Aft  has  essentially  achieved 
the  value  associated  with  a d i ffract ion- I imi ted  transfer  function,  we 
intuitively  expect  the  s igna I -to-no i se  ratio  Q to  increase  in  propor- 
tion to  l7  for  higher  light  levels,  in  accord  with  the  Poisson  character 
of  the  noise.  Unfortunately,  when  a flat-topped  S is  used,  Q ulti- 
mately increases  as  , a consequence  of  the  fact  that  the  restoration 

filler  places  U<o  great  an  emphasis  un  restoring  frequency  components  that 
are  extremely  close  to  the  cutoff  frequency  To  complicate  matters 

even  further,  our  previous  plots  of  Q were  found  to  ultimately  increase 
as  hP  , rather  than  N . This  behavior  was  found  to  arise  from  our 
use  of  an  approximate  expression  for  the  diffraction  limited  OTF  [1,  Eq . ( 59) J 
which  approaches  zero  at  a different  rate  than  the  exact  expression.  To 
remove  this  pathological  dependence  of  Q on  the  behavior  of  the  OTF 
infinitesimally  close  to  cutoff,  we  have  found  it  highly  desirable  to 

A 

use  the  S of  Eq.(5)-  This  choice  has  the  added  advantage  that  the 
Ideally  filtered  object  will  always  be  non-negative,  whereas  this  is  not 
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necessarily  the  case  with  a flat-topped  S. 

The  blur  transfer  function  assumed  here  (as  in  the  previous  report) 
is  of  the  form 

B(fl)  - BT(S3)BA(Q)  (6) 

where  B^.(^)  is  the  optical  transfer  function  of  the  telescope,  while 
BA(ft)  is  the  optical  transfer  function  of  the  atmosphere.  For  simplicity, 
we  assume  that  the  telescope  has  a clear  circular  aperture  (no  central 
obscuration)  and  is  diffraction  limited,  in  which  case 


(7) 


for  <_  Mg.  For  the  atmospheric  OTF,  we  assume  that  the  exposure  time 
is  long  compared  to  the  atmospheric  fluctuation  time,  and  use  an  OTF  of 
the  form 


for  il  <_  Here  D is  the  diameter  of  the  telescope  aperture  and  a 

takes  on  the  values  0,  H or  l as  follows: 
a * 0:  No  t i I t removal ; 

a ■ H'<  Perfect  tilt  removal,  far  field  atmospheric 
propagat ion; 

a ■ I:  Perfect  tilt  removal,  near  field  atmospheric 
propagat ion. 

Tilt  removal  refers  to  the  case  of  a telescope  with  a servo-controlled 
mirror  which  removes  the  linear  terms  In  a least-squares  fit  to  the 
atmospherically  distorted  wavefront.  The  result  Is  nearly  (but  not 
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exactly)  the  same  as  using  a perfect  centroid  tracker  to  determine  the 
requ i red  mi rr or  tilt. 

In  Fig.  2a,  b and  c,  we  present  the  plots  of  restored  bandwidth 

Afi  vs.  N for  the  three  cases  described  above.  The  curve  for  r ■ » 

o 

represents  the  ideal  behavior  of  a diffraction-limited  telescope  in  the 
absence  of  atmospheric  turbulence.  In  all  cases  the  telescope  mirror 
diameter  was  assumed  to  be  1.524  meters  (60  inches)  and  the  mean  wave- 
length was  taken  to  be  5*10  ^ meters.  For  large  enough  N , all  of  the 
curves  eventually  reach  the  d i f f ract ion- 1 imi ted  performance  curve 
(r  = »),  but  for  the  range  of  N shown,  only  the  rQ  = 20cm  curve  of 
Fig.  2(c)  accomplishes  this. 

Also  shown  on  each  drawing  are  three  dotted  horizontal  lines.  These 
lines  represent  the  value  of  AQ  for  rQ  ■ 5cm,  10cm  and  20cm  when  n£ 
post  processing  is  applied  to  the  detected  image  (i.e.,  H - 1 for 
Q <_  ftp  , zero  otherwise).  By  comparing  the  corresponding  dotted  and 
solid  lines,  it  is  possible  to  deduce  quantitatively  just  how  much  image 
bandwidth  has  been  gained  by  using  the  least-squares  restoration  filter. 

Figures  b and  c present  plots  of  the  signal-to-noise  ratio  Q 
vs.  N"  , again  for  four  values  of  rQ . eventually,  all  of  these  curves 
bend  upwards  and  increase  in  proportion  to  N , running  parallel  to,  but 
below,  the  rQ  ■ 00  curve.  However,  for  the  range  of  N shown,  this  upward 
turn  can  be  seen  only  for  the  rQ  ■ 20cm  curve  of  Fig.  3(c).  By  comparing 
Figs.  2(c)  and  3(c),  we  note  that  the  upward  turn  in  Q takes  place  at 
values  of  N where  the  restored  bandwidth  is  closely  approaching  the 
d 1 ff  ract  ion- 1 imi  ted  curve  (rQ  • “>) . Thus,  once  the  restoration  filter 
has  succeeded  in  restoring  frequency  components  to  the  level  appropriate 


for  a diffraction-limited  telescope,  ail  further  increases  of  N have 
the  effect  only  of  increasing  the  s ignal -to-photon-no i se  ratio,  with  a 
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resulting  increase  of  Q in  proportion  to  N.  For  the  rQ  - 10cm, 

perfect  tilt  removed,  near  field  curve,  this  transition  has  been  found 

— 21  — 
to  take  place  at  N " 10  . which  is  far  off  the  scale  of  N shown. 

i V.  PERFORMANCE  OF  A SYSTEM  WITH  BOTH  PREDETECTION 
COMPENSATION  AND  PQSTDETECT I ON  COMPENSATION 

We  consider  now  the  photon-limited  performance  of  a system  which 
utilizes  both  predetection  compensation  and  postdetection  compensation. 
With  reference  to  Fig.  A,  a portion  of  the  incoming  light  is  sent  to  a 
wavefront  sensing  system,  which  measures  the  atmospheric  wavefront  dis- 
tortions in  real  time.  The  wavefront  sensor  then  drives  a wavefront 
corrector  (usually  a deformable  mirror),  which  removes  the  atmospheric 
distortions  from  the  incoming  light.  In  addition,  a portion  of  the  light 
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Fig.  4.  Compensated  imaging  system  block  diagram 
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leaving  the  wavefront  corrector  is  sent  directly  to  an  image  detector. 

The  detected  image  is  then  subjected  to  post-processing,  following  which 
the  final  image  is  displayed. 

Two  types  of  error  arise  in  this  system  as  a result  of  the  finite 
photon  flux  available  at  the  telescope.  First,  the  wavefront  sensor 
functions  imperfectly  in  the  presence  of  photon  noise,  with  the  result 
that  to  some  extent  the  wavefront  is  improperly  corrected,  and  the 
recorded  image  suffers  residual  blur.  Secondly,  the  finite  photon  flux 
sent  to  the  image  detector  limits  the  ability  of  the  post-processor  to 
remove  these  residual  defects. 

Our  goal  is  to  study  image  qualt  /,  as  fined  by  the  two  para- 
meters Afi  and  Q of  the  previous  se<  mo  , as  a function  of  the  mean 
number  of  photoevents  N supplied  by  ‘ lelescope.  The  post-processor 
will  again  be  a linear  invariant  lep  squares  filter.  In  fact,  the 
entire  formalism  represented  by  Fig.  1 can  again  L sed,  but  with  a blur 
transfer  function  appropriate  for  the  partially  precompensated  system. 

Our  model  for  the  partially  compensated  wavefronts  is  that  of  a 
spatially  stationary  gaussian  random  phase  screen.  This  model  is,  of 
course,  only  an  approximation,  for  the  fixed  locations  of  the  discrete 
actuators  on  the  deformable  mirror  introduce  a spatial  non-stat ionarity, 
sometimes  referred  to  as  the  "waffle  effect".  However,  the  statistic- 
ally stationary  phase  screen  provides  a reasonable  first  order  model. 

From  such  a model,  It  is  easily  shown  that  the  blur  transfer  function 
associated  with  the  partially  compensated  wavefronts  is  of  the  form 

Bc(n)  - expj-  i D^(Ah)j  (9) 

where  X is  the  mean  wavelength  and  D (r)  is  the  structure  function  of 

4> 


the  residual  phase  errors, 


Vr)  " E[(  4>2)2]  (10) 

Here  4> ^ and  are  two  values  of  the  phase  at  points  separated  by 

distance  r , and  it  has  been  tacitly  assumed  that  the  structure 
function  ic  circularly  symmetric.  Implicit  in  this  model  is  the  assump- 
tion of  "near  field"  atmospheric  propagation,  i.e.,  intensity  scinti- 
llation effects  are  negligible. 

In  evaluating  the  performance  of  the  overall  system,  we  have  found 
it  necessary  to  specify  the  character  of  the  precompensation  system  more 
precisely.  We  assume  here  that  che  wavefront  sensor  consists  of  a 
shearing  interferometer  that  produces  an  array  of  noisy  phase-difference 
measurements.  These  measurements  are  applied  to  a phase  reconstruction 
network,  which  produces  an  array  of  noisy  phase  estimates  from  the  phase 
diffe-ence  measurements.  The  phase  estimates  are  then  used  to  drive  a 
deformable  mirror.  In  general,  the  phase  estimates  are  based  on  a measure- 
ment interval  that  is  some  small  fraction  of  the  total  Integration  time 
used  to  detect  the  corrected  image.  In  the  specific  results  presented 
later  we  have  assumed  a ratio  of  total  image  integration  time  T to 
phase  measurement  time  x of  loV  This  ratio  is  appropriate,  for  example, 
when  T ■ 1 second  and  x ■ 100  microseconds.  We  shall  further  assume 
that  any  errors  Induced  by  the  fact  that  the  active  mirror  is  driven  to 
correct  for  errors  that  were  measured  a short  time  in  the  past  are  ncgllgiole. 
In  other  words,  we  assume  that  the  delay  between  measurement  and  correction 
is  negligibly  small  compared  to  the  other  errors. 
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The  phase  structure  function  consists  of  two  parts, 


Mr)  - Djf)(r)  + Djn)(r) 
9 9 9 


00 


,(0 


where  D^’'(r)  is  the  structure  function  component  due  to  fitting  error 
and  (r)  is  the  structure  function  component  due  to  noise.  McGlamery 

[2]  has  estimated  these  structure  functions  by  simulation  of  a system 
having  317  subapertures  within  a 1.6  meter  aperture.  In  turn,  we  have 
fit  McGlamery's  simulation  results  with  analytic  expressions.  tye  find 
excellent  fit  to  the  simulation  data  by  the  following  expressions: 


D^(Ih)  = 2o* 


D^n)(\n)  « 2a* 

9 n L 


/(*•*  kl 
-hk) 
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2 2 

Here  and  an  are  the  phase  variances  due  to  fitting  error  and  noise, 

respect i vel y. 

2 

Using  Hudgin's  [3]  results,  we  can  express  as 


or  - 0.31 6(  — 

f Vo 
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(rad)' 
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where  L is  the  spacing  of  the  actuators.  With  K subapertures  in 

Q 

an  unobstructed  aperture  of  diameter  D meters,  we  have  approximately 
that 


of  - 0.316 


^5/3 


,/K1 


(14) 


Turning  attention  to  the  phase  variance  due  to  noise,  it  is  well 
established  that  the  variance  of  the  electronic  phases  measured  by  the 
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wavefront  sensor  are 
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where  y is  the  electronic  fringe  contrast  and  Ng  is  the  average 

number  of  photoevents  per  phase-difference  sensor.  For  a point-source 

2 

object,  we  assume  y ■ 1.  Taking  account  of  the  fact  that  there  are 
two  shears  involved,  one  with  respect  to  x and  the  second  with  respect 
to  y , we  have 


(16) 


where  g is  the  fraction  of  photoevents  sent  to  the  wavefront  sensor 
(0  < g < I),  and  N is  the  total  number  of  photoever.ts  collected  by  the 


telescope.  (Implicit  in  our  expressions  are  the  assumptions  of  equal 
losses  and  equal  quantum  efficiencies  in  the  wavefront  sensing  and  image 
detecting  branches  of  the  optical  system.)  Thus,  for  a point  source 
object, 


M (rad)2 

gNx 


(17) 


The  variance  of  the  phase  difference  measurements  is  related  to 


through 


% • 


(18) 


where  S Is  the  shear  distance  and  Lg  is  the  sensor  spacing.  For  a 

point  source  object,  we  assume  that  the  shear  distance  can  be  made  as 

2 2 

big  as  the  sensor  spacing,  yielding  ■ ae. 

We  use  Fried's  results  [4]  to  relate  the  variance  of  the 

estimated  phases  to  the  variances  of  the  phase  differences, 


o2  - 0.655811  + 0.2444  £n  K]  a1..  . (19) 

n ay 
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Combining  (19).  (18),  and  (17),  we  find 


2 

on 


0 . 6558 [ 1 + 0 . 2444  in  K] 


(20) 


Equations  (9),  (12),  (14)  and  (20)  now  provide  us  with  a complete 
specification  of  the  optical  transfer  function  of  the  precompensation 
system.  This  transfer  function  must  be  multiplied  by  the  telescope 
transfer  function  By(£2)  to  obtain  the  entire  blur  transfer  function. 

Finally,  we  take  into  account  the  effect  of  the  post-compensation 
system.  Noting  that  the  average  number  of  photoevents  appearing  in  the 

— A A 

image  is  ( 1 - 0 ) N , and  using  S = By  for  the  ideal  filter,  we  use 
Eq. (47)  of  Ref.  1 to  write  the  product  of  blur  and  deblur  OTF's  as 


[1  + ( 1 -B) N) By  B* 

BH  = , (21) 

1 + (l-B)N  By  B^ 

A 

where  a point-source  object  (<£Q  » 1)  has  been  explicitly  assumed. 
From  this  quantity,  the  image  quality  parameters  A£2  and  Q can  be 
calculated  by  numerical  integration, 


and 


Aft 


n|B(n)H(n)  dn 


(22) 


00  fi  a rt 

C By(n)B*(n)  Qdfi 

_ l i + (i-e)N  b?({j)bJ(o) 

Q - (l-B)Njj 1 5 • (23) 

j*  By  (fi)  Q d St 

1 i + 0-e)N  B*(n)B*(n) 

We  now  present  results  of  computer  calculations  of  AH  and  Q 
based  on  the  theory  outlined  above.  In  all  cases  a single  point-source 
object  Is  assumed,  and  the  number  of  wavefront-sensor  subapertures  is 
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I Ocm 


317.  Figures  5a  and  b show  plots  of  A?5  and  Q vs.  N for  rQ  - 

and  three  different  values  of  B (0.3,  0.9,  0.99).  Note  that  the 

_ 6 7 

diffraction  limited  bandwidth  is  reached  for  N in  the  range  10  to  10  , 
depending  on  3.  The  values  of  Q achieved  are  highest  for  large  3 
when  N is  below  the  value  which  achieves  diffraction  limited  perform- 
ance. However,  for  values  of  N above  this  level,  small  values  of  3 
are  superior.  Once  diffraction  limited  bandwidth  has  been  reached, 
there  is  no  point  to  sending  more  light  to  the  wavefront  sensor.  Hence 
smaller  and  smaller  fractions  of  the  incident  light  should  be  sent  to 
the  wavefront  sensor  as  N increases  beyond  its  critical  value. 

The  optimization  of  the  splitting  ratio  3 is  explored  in  more 
detail  in  Figs.  6a  and  b.  Aft  and  Q are  shown  as  functions  of  3 
for  various  values  of  N , with  rQ  fixed  at  10cm.  The  optimum  value 
of  3 is  best  determined  from  Q,  in  Fig.  6b,  for  Aft!  ultimately 
saturates  at  the  d i ff ract ion- l imi ted  bandwidth,  while  Q has  no  such 
saturation  limit.  The  optimum  vaiue  of  3 is  seen  to  depend  rather 
strongly  on  the  object  brightness  N.  For  hT  small  (<  10^)  , the 
optimum  3 is  close  to  unity,  indicating  that  most  of  the  light  should 
be  sent  to  the  wavefront  sensor.  However,  for  N greater  than  about 
IQ7,  the  optimum  value  of  3 shifts  abruptly  to  rather  small  values, 
Indicating  that  most  of  the  light  should  be  sent  to  the  image  detector. 
This  behavior  simply  indicates  once  again  that  when  A?2  has  reached 
the  diffraction  limit,  further  increases  of  light  flux  should  be 
allocated  to  the  Image  detector  for  best  performance. 

Figures  7a  and  b show  the  dependence  of  system  performance  on 
the  atmospheric  coherence  diameter  r . In  this  case  3 Is  fixed  at 
0.9;  Afi  and  Q are  plotted  vs.  N.  Concentrating  on  Fig.  7b,  for  F 
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ft 


Fig.  6.  Restored  bandwidth  (a)  and  quality  factor 
(b)  vs.  3 for  ro-IOcm  and  various  values 
of  N. 
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below  the  critical  range  of  values  that  produce  d i f f ract ion-1 im i ted 
bandwidth,  the  difference  of  values  of  Q achieved  as  rQ  varies  is 
small.  For  N above  this  critical  range,  the  change  in  Q as  rQ 
varies  is  larger,  but  less  than  an  order  of  magnitude. 

In  the  final  figure  (Fig.  8),  we  investigate  the  effects  of 
removing  certain  parts  of  the  overall  pre-  and  post-compensated  imaging 
system.  Plots  of  Aft  vs.  N are  shown,  with  rQ  fixed  at  10cm. 

Curves  are  shown  for  pre-  and  post-compensation  (0  = 0.9),  pre-compen- 
sation alone  (0  » 0.9),  post-compensation  alone  (0  *1.0,  with  and 
without  simple  tilt  removal).  A number  of  observations  can  be  made  from 
these  curves.  First,  comparing  the  curve  for  pre-compensation  only  (b) 
with  the  curve  for  both  pre-  and  post-compensation  (a),  we  see  the  differ- 
ences are  greatest  in  the  transition  region  where  the  restored  bandwidth 


4 e • io 

Ioqiq  R 


Fig.  0.  Comparison  of  AH  vs.  N*  curves  for 
various  compensation  systems. 

(a)  Precompensation  and  Postcompensation. 

(b)  Precompensat Ion  Only . 

(c)  Postcompensation  and  Tilt  Removal, 
Near  Field. 

(d)  Postcompensation,  Long  Exposure. 

(e)  Tilt  Removal  Only,  Near  Field. 

(f)  l(o  Compensation,  Long  Exposure. 
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is  approaching  saturation.  For  larger  values  of  N , the  difference 

between  the  two  saturation  levels  is  attributable  to  fitting  errors.  In 

one  case  the  effects  of  fitting  errors  are  post-compensated,  in  the  ^ 

second  case  they  are  not. 

Next  we  note  that  over  a large  range  of  N , the  values  of  Afi 
achieved  with  post-compensation  only  fall  far  below  those  achieved  with 
pre-compensation  only.  This  conclusion  remains  valid  whether  overall 
tilt  removal  is  used  or  not  in  the  post-compensation  case.  Note,  howe'er, 
that  if  N is  sufficiently  small,  ultimately  the  pre-compensation 
system  begins  to  degrade  the  performance,  for  the  wavefront  measurements 
become  so  noisy  that  they  add  errors  that  are  worse  than  those  intro- 
duced by  the  atmosphere.  Thus,  for  example,  with  N « 10**  , best  per- 
formance is  achieved  by  turning  the  pre-compensation  system  off  and 
performing  only  post-compensation.  However,  the  values  of  AS?  achieved 
in  this  regime  are  not  much  better  than  what  could  be  achieved  with  the 
unmodified  original  telescope  without  any  post-compensating. 

V.  SENSITIVITY  OF  SPECKLE  IMAGING 

By  the  term  "speckle  imaging"  we  mean  the  technique  Introduced  by 
Labeyrie  [5]  which  depends  on  averaging  the  spectral  density  distributions 
of  a sequence  of  k short-exposure  photographs.  If  Ij (v)  represents 
the  Fourier  transform  of  the  short-exposure  image  lj(x)  , normal- 

ized to  have  value  unity  at  the  origin,  then 

Ij(v)  “ Bj(v)0(v) 

where  Bj  (v)  represents  the  optical  transfer  function  (OTF)  of  the  k 

a 

exposure,  and  0(v)  represents  the  normalized  Fourier  spectrum  of  the 
object  (assumed  identical  for  all  k measurements). 
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The  Labeyrie  technique  uses  < independent  measurements  of 
, |I  (v) I to  attempt  to  estimate  |0(v)|  . The  basis  for  this  method 

i rests  upon  the  relation 

^|ij(v)  |2^  - <^|  Bj  (v)  | ^ j 0(v)  | 2 (24) 

and  the  fact  that  [6],  in  the  mid-frequency  region  of  the  diffraction- 
1 imi ted  OTF, 

(|'jmi2)s(t)  btm  (25) 

where  rQ  is  Fried's  atmospheric  coherence  diameter;  D is  the 
diameter  of  the  telescope  used;  BT(v)  represents  the  diffraction- 
limited  OTF  of  the  telescope;  and  it  has  been  assumed  that  D/rQ  >>  1. 
Thus  in  the  mid-frequency  range  of  the  telescope  passband, 

r 2 

<^|ij(v)|2^-  •£  Bt(v)  |0(v)  | 2 (26) 

i where  D , rQ  and  B^.  are  presumed  known,  |I(v)|  is  to  be  estimated 

' / * 2v 

from  the  sequence  of  < pictures,  and  (|0(v)|  / is  the  ultimate  quantity 
of  interest. 

For  any  particular  sequence  of  k pictures,  the  finite  average 

I < 2 / 2\ 

— I | I . (v) | is  only  approximately  equal  to  the  true  average  (|I.(v)|  > 
j K J-l  J \ J / 

The  statistical  fluctuations  of  the  finite  average  arise  from  two  sources: 

(1)  The  statistical  fluctuations  of  Bj (v)  from  picture  to 
picture,  as  caused  by  the  changing  atmospheric  turbulence, 

and 

(2)  The  statistical  fluctuations  of  ij(x)  caused  by  the  finite 
number  of  photcavents  detected  In  each  picture. 
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Both  of  these  sources  of  fluctuations  lead  to  errors  in  our  estimate  of 


<!¥'■)!> 

The  effects  of  these  two  types  of  fluctuations  on  an  estimate  of 
<(|i<v)l2)  have  been  evaluated  in  our  previous  technical  report  [1]. 

From  Eq.(3l)  of  that  report,  we  have  that  the  r .m. s s ignai -to-noi se  ratio 
associated  with  the  measurement  is,  for  a single  frame  (<  * 1 ) , 


(v)  +-~<fij(\>)  + 

N 


♦i 

(v) 

(N)2 

^1  + 4$.  (v)  + . (2v) 

* -M- 

(N)3/ 

(27) 

- 


where  N Is  the  average  total  number  of  photoevents  detected  in  a 
single  frame,  and 

41  i (v)  ■ ^|I(v)|2^  • (28) 

A similar  expression  has  been  derived  by  Roddier  [7].  When  the  measure- 
ment is  based  on  k independent  frames  of  data,  the  S/N  is  increased 
by  Jp  , 


The  frequency  range  of  greatest  interest  is  at  the  high  end  of  the 
diffraction-l imi ted  telescope  passband.  Hence,  If  vQ  represents  the 
d i f f ract ion- 1 imi ted  cutoff  frequency  of  the  telescope,  we  assume  that 
v > vQ/2  , in  which  case  the  term  d> j (2v)  of  Eq. (27)  can  be  replaced  by 
zero.  For  simplicity  of  notation  we  define 

i0(v)  - | 0(v) | 2 (30) 

to  represent  the  normalized  spectral  density  of  the  object  itself,  in 
which  case  Eq.(26)  becomes 
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(v) 


(v)<f>Q  (v) 


Finally,  we  define  the  quantity 


fr)! 


(3D 


(32) 


which  represents  the  average  number  of  photoevents  per  atmospheric 
coherence  area.  This  number  can  equivalently  be  regarded  as  the  average 
number  of  photoevents  per  speckle  in  the  image. 

With  the  above  definitions  and  some  simple  algebraic  manipulations, 
we  obtain  from  Eq.(27)> 


{['  + n 8t 


(v)<j>Q(v) 


T~T  “ " TV 

— 1 + ^n  BT (v) 4?  (v)  l 

N L JJ 


(33) 


valid  for  v > v0/2. 

Three  important  regions  of  this  relationship  can  be  identified,  in 
each  of  which  a different  dependence  on  n is  found 

(1)  For  n BT<f>Q  >>  1 , 


1 (independent  of  n)  ; 

(2)  n Bt*q  <<  1 but  N >>  1 , 

(n)  * n Bj(v)^0(v)  ; 

A A 

(3)  n By(J>0  <<  1 but  N <<  I , 

(I)0’  ' (n)V2  (^)  iT(v)J0(v)  . 

t 

Two  important  points  should  be  emphasized.  First,  the  maximum  possible 
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single-frame  signal -to-nol se  ratio  for  speckle  i nterfero.  ietry  is  unity, 
regardless  of  how  bright  the  object  may  be.  Second,  only  when  the  average 
, number  of  phof.oevents  per  picture  is  much  less  than  unity  does  increasing 

i the  telescope  aperture  diameter  improve  the  s i gna l -to-no i se  ratio. 

The  single-frame  s i gna 1 -to-noi se  ratio  is  plotted  in  Fig.  9 for 
the  special  case  D ■ 1.52  meters,  rQ  » 10cm,  v/vQ  = 0.8  and  ■>  1 . 


Fig.  9*  Single-frame  rms  s ignal -to-noi se  ratio 
vs.  average  number  of  photoevents  per 
r0-cell,  speckle  Imaging. 

The  thr^e  separate  regions  of  dependence  on  IT  are  clearly  distinguish- 
able. 
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METRIC  SYSTEM 


BASE  UNITS: 


^jant'ty 

_ Unit 

SI  Symbol 

Formula 

length 

metre 

m 

mas* 

kilogram 

kg 

time 

second 

s 

electric  currant 

ampere 

A 

thermodynamic  temperature 

kelvin 

K 

amount  of  aubetance 

mole 

mol 

lumlnoua  intenaity 

candela 

cd 

SUPPLEMENTARY  UNITS; 

plane  angle 

radian 

rad 

aolld  angle 

steradlan 

•r 

... 

DERIVED  UNITS: 

Acceleration 

metre  per  second  squared 

m/a 

activity  (of  a radioactive  aource) 

disintegration  per  second 

(disintegration)/! 

angular  acceleration 

radian  per  second  squared 

rad/s 

angular  velocity 

radian  per  second 

rad/s 

area 

square  metre 

m 

denalty 

kilogram  per  cubic  metre 

kg/m 

electric  capacitance 

farad 

F 

A-s/V 

electrical  conductance 

siemens 

S 

W 

electric  field  atrength 

volt  per  metre 

v.'m 

electric  inductance 

henry 

H 

V./A 

electric  potr  .dial  difference 

volt 

V 

W/.A 

electric  rnaiatance 

ohm 

V 

electromotive  force 

volt 

V 

W/A 

energy 

joule 

) 

N-m 

entropy 

joule  per  kelvin 

J/K 

force 

newton 

N 

kg-m/s 

frequency 

hertz 

Hz 

(cycle)/! 

Illuminance 

lux 

lx 

Im/m 

luminance 

candela  per  square  metre 

cd/m 

luminoua  flux 

lumen 

!m 

cd-ar 

magnetic  field  stren  a 

ampere  per  metre 

* 

A/m 

magnetic  flux 

weber 

Wb 

V-a 

magnetic  flux  density 

tesla 

T 

Wb/m 

magnetomotive  force 

ampere 

A 

power 

watt 

W 

j/a 

pressure 

pascal 

Ps 

N/m 

quantity  of  electricity 

coulomb 

C 

A-s 

quantity  of  heat 

joule 

1 

N-m 

radiant  intensity 

watt  per  steradlan 

W/sr 

specific  heat 

joule  per  kllogram-kelvln 

|/kg-K 

stress 

pascal 

Pa 

N/m 

thermal  conductivity 

watt  per  matre-kelvin 

W/m-K 

velocity 

matre  per  second 

m/s 

viscosity,  dynamic 

pascal-second 

Pa-s 

viscosity,  kinematic 

square  metre  per  second 

m/a 

voltage 

volt 

V 

W/A 

volume 

cubic  metre 

m 

we.  -number 

reciprocal  metre 

(wave)/m 

work 

joule 

i 

N-m 

SI  PREFIX  SS: 


Multiplication  Factor! 


Prefix  SI  Symbol 


1 000  000  000  000  = to" 

1 000  000  000  - 1 O’ 

1 000  000  - 10* 

1 000  - 10’ 
100  - 10' 

10  » 10' 
0.1  *=  10-' 
0.01  = 10-* 
0.001  - 10-' 
O.tKJO  001  » 10- * 
0.000  000  001  =■  10- * 
0.000  000  000  001  ' 10* '« 
0.000  000  000  (too  001  - 10-" 

0.000  000  OOO  000  000  001  10"  •» 

To  ba  avoided  whera  poesible 


lore 

■Ida 

mega 

kilo 

hacto* 

daka* 

decl* 
cent  I" 
mllll 

micro 

nano 

filoo 
amlo 
alto 
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